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ABSTRACT: The paper presents the results from an analysis of the correlation and 
regression, and assessment of the recurrence of correlation of 13 parameters determined in the 
soils of five regional varieties of Tilio-Carpinetum forests. It has been found that: (1) most of 
the parameters considered are positively, linearly correlated; (2) the interrelationships of most 
of the characters analysed provide a better descritpion of the soil environment of the particular 
regional varieties of Tilio-Carpinetum forests than of the habitat type analysed. 
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1. INTRODUCTION 


In the relevant literature, especially that of the last two decades, much attention has 
been given to the stochastic relationships of the parameters determining the 
physicochemical properties of the soil (Canarache 1968(1967), Banin and 
Amiel 1969(1970, Fassbender etal. 1969(1970), Singh ans Shah 
1979, Chudecki and Błaszczyk 198, Czarnowska 1980). 
Studies of this type were done while carrying out standard analyses of soil properties, as 
well as in order to find coefficients to be used in the calculation of other parameters, 
sometimes difficult to determine in an empirical way (Canarache et al. 1969 
(1968, Eh walTd and Krüger 1976 Gebhard and Illner 1976, 
Trzecki 1976, Turski etal. 1978, Urvas etal. 1978 Mckeaque 
and Wolynetz 1979, Gawlik and Zawadzki 1980, Glinski 
and Stepniewska 1980) 


[283] 


284 Urszula Wierzchowska 


The occurrence of stochastic relationships among many soil parameters and the 
types of these associations have been fairly well documented in the literature. 
Nevertheless, relationships found for two characters on one occasion are not always 
confirmed by results from investigations carried out for other soils, sometimes with 
similar physicochemical properties. For this reason, in the correlation and regression 
analyses attention was paid not only to the degree and type of a relationship between 
two variables, but also to their recurrence. 

The aim of the present study was to carry out a statistical analysis of the 
relationships between selected, ecologically significant soil properties exemplified by 
the soils of mixed oak-hornbeam forests of the Tilio-Carpinetum type. This study is 
part of more comprehensive investigations into the soil-condition diversity of these 
communities (Wierzchowska 1981). 

Mixed oak-hornbeam communities of the Tilio-Carpinetum type represent 
mesotrophic, multi-species deciduous forests of the eastern part of Poland. A high 
floristic diversity of the oak-hornbeam phytocoenoses and its clearly regional nature 
are “reflected” by the typological differentiation of this unit into regional varieties, 
namely: the sub-Boreal variety, the central-Polish variety, the Little-Polish variety with 
two altitudinal forms, the Volhynian and the Nida varieties (Matuszkiewicz 
and Matuszkiewicz 1981). 

The soil environmental conditions of mixed oak-hornbeam forests of the Tilio- 
Carpinetum type are characteristically highly variable. Within the ecological ampli- 
tude of these communities statistically significant regional soil condition variation has 
been found. The biotope types distinguished correspond with the above-enumerated 
typological units of oak-hornbeam forests (Wierzchowska 1981). 

Data relating to different biotope types, representing the same habitat type, have 
been used in correlation and regression analyses of some soil parameters. These 
characters are most often taken into account in ecological studies when the soil 
environment is analysed as one of the most important factors determining vegetation 
diversity. 

The analysis of recurrence of soil parameter correlation in the five biotope types 
taken into account is a trial for answering the question concerning the usefulness of the 
results from correlation and regression analyses in the description of the soil 
environment of plant communities. 


2. MATERIAL AND METHODS 


The basic material of the paper is numerical data on 13 characters aescribing tne 
properties of the humus horizon of the soils of Tilio-Carpinetum phytocoenoses. A 
total of 492 soil samples have been analysed. The samples were collected during the 
growing seasons of 1977 and 1978 in the whole area of occurrence of these communities 
in Poland. 

The data taken into account include: per cent content of particles below 0.02 mm, 
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amount of available phosphorus in mg per 100 g of soil, content of metal cations 
(Ca**, Mg**, K*, Na*) in me: 100 g soil ', hydrolytic acidity (H,) in me: 100 g 
soil *, pH in H,O, pH in KCI, per cent content of organic carbon (C org.) and metal 
cation total (S) in me- 100 g soil” *, ion exchange capacity (T) in me: 100 g soil * and 
degree of saturation of ion exchange capacity with metal cations (V.). These characters 
have been determined by using the methods generally accepted in soil science 
(Swiecicki etal. 1966). 

For the purpose of statistical analysis, all the above parameters were by assumption 
considered equivalent characteristics of the soil environment. 

The degree of correlation between particular soil parameters is described by the 
linear correlation coefficient (r) and the parabolic correlation coefficient (R), calculated 
for all the character pairs analysed. 

Linear regression equation y = ax + b and parabolic regression equation 
y = ax? + bx + c describe the type of relationships between the soil properties studied. 
Values of the a, b and a, b, c parameters have been estimated by the least square method. 
The regression significance hypothesis has been verified by comparing empirical r and 
R values with critical values from tables (Elandt 1964). Statistical inferring was 
done at the 5%, error risk. 

Correlation and regression analyses were carried out separately for numerical data 
characterizing the biotopes of five mixed oak-hornbeam forest units, namely: the 
sub-Boreal variety (I), central-Polish v. (II), upland form of the Little-Polish v. (IID), 
piedmont form ofthe Little-Polish v. (IV) and the Volhynian variety (V). For each of the 
five combinations considered the degree of correlation of 78 character pairs has been 
analysed. 

On the basis of the number of statistically significant correlations, found for a 
specific character pair in all combinations analysed, a diagram has been drawn to show 
the degree of correlation of the properties of the soils studied. Lines drawn on the basis 
of regression equations represent the interrelationships of the characters studied. 


3. RESULTS AND DISCUSSION 


The results of the correlation analysis indicate that most of the soil parameters 
taken into account are significantly correlated. The list of associated soil character 
pairs and the number of these pairs are more or less different in each of the oak- 
hornbeam varieties (Tables 1 — 5). For the soils of the sub-Boreal variety 46 correlated 
character pairs have been found, for oak-hornbeam forests of the central- Polish variety 
— 51 (therein 2 character pairs specific only to it), for oak-hornbeam forests of the 
upland form of the Little-Polish variety — 41 (therein | pair specific only to it), for the 
piedmont form of this variety — 40 (therein 4 pairs typical of it), and for the soil of 
oak-hornbeam forests of the Volhynian variety — 55 pairs of correlated characters 
(therein 2 pairs of characters specific only to it). 
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Table 2. Values of coefficients of linear (r) and parabolic (R) correlation of soil properties of the central-Polish variety of Tilio-Carpinetum (limit value: 
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Table 3. Values of coefficients of linear (r) and parabolic (R) correlation of soil properties of the upland form of the Little-Polish variety of Tilio-Carpinetum 
(limit value: rg gs = 0.325, Ry ys = 0.400) 
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Table 4. Values of coefficients of linear (r) and parabolic (R) correlation of soil properties of the piedmont form of the Little-Polish variety of Tilio-Carpinetum 
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Table 5. Values of coefficients of linear (r) and parabolic (R) correlation of soil properties of the Volhynian variety of Tilio-Carpinetum (limit value: 
Pu ns = 0.444, Ry gs = 0.545) 
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Fig. 1. Diagram presenting statistically significant correlations of the Tilio-Carpinetum forest soil properties 
studied 
m.f. — mechanical fraction, a — sub-Boreal variety, b — central-Polish variety, c — upland form of the 
Little-Polish variety, d — piedmont form of the Little-Polish variety,e — Volhynian variety 


Significant correlation, of the parabolic type, in all the five biotope types has been 
found for 27 correlative pairs, in four biotope types — for 10 pairs, in three types — for 
9,in two — for 10, in one biotope type — for 11 pairs, and a lack of correlation has been 
found for 11 character pairs. 

In the case of linear correlation in all the five biotope types a significant correlation 
has been found for 27 correlative pairs, in four biotope types — for 11, in three types — 
for 5,in two — for 13, in one biotope type — for 12 pairs, and a lack of correlation has 
been found for 10 character pairs (Fig. 1). 

The number of correlated pairs common to the particular oak-hornbeam varieties 
(I-II, I—HI, I—IV, I— V, II— III, I—IV, W—V, HI—IV, III— V, IV—V) varies 
considerably. The largest number, as many as 47, of pairs has been found for the 
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central-Polish and the Volhynian varieties (II — V), and the smallest number — 30 — 
for the central-Polish variety and the piedmont form of the Little-Polish variety (11 
— IV). Apart from this, for the first variety pair the largest number, as many as 7, of pairs 
of correlated characters have been found, for which statistically significant correlation 
has only been found in the case of this pair of oak-hornbeam varieties. The number of 
character pairs common to the remainder of the variety pairs ranges from 33 to 40, 
whereas the number of character pairs typical of only one specific variety pair is 2, at the 
most. 

The number of correlations (linear and parabolic) found in the five combinations 
analysed has been adopted as a measure of correlation of two characters (Fig. 2). 
Characters for which statistically significant correlation has been found in at least eight 
(out of ten possible) cases have been considered highly correlated. The adoption of such 
an assumption makes it possible to highlight those parameters for which recurrence of 
a statistically significant interrelationship has been found. 

A significant correlation of two characters found only for one combination may 
characterize an association of soil parameters in the given system of abiotic conditions, 
or the sample alone. 
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Fig. 2. Diagram showing the 
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On the basis of the analysis carried out the following have been classified as highly 
correlated: metal cations, total cation content, degree of ion exchange capacity 
saturation and the sorptive capacity of the soil. The latter is further correlated with the 
content of organic matter expressed by the amount of organic carbon. Metal cations, 
except potassium, and the degree of ion exchange capacity saturation with cations are 
also correlated with the soil pH. For potassium a close correlation with silt aee clay 
content has been found. Of the cations analysed most highly correlated with the 
content of organic carbon was sodium. 

A character showing no correlation with other soil properties is the content of 
available phosphorus. A poor correlation with the parameters studied has also been 
found for the hydrolytic acidity. It shows the highest, relative, correlation with the 
content of organic carbon and with the soil sorptive capacity. 

The association of characters found in the present study is confirmed by the papers 
ofotherauthors (Markert 1965, Russell and More 1967, Cloos et 
al. 1968, Hudcova and Kovarova 1969, Matuszkiewicz 1974, 
Pavlov and Kjuleva 1976). 

The groups of highly correlated characters that can be distinguished in the soil 
parameter correlation diagram are associated with the properties of the material 
constituting the parent rock of the soil, or directly with the features of the soil 
mineral-organic complex (Russell and More 1967, Matuszkie- 
wicz 1974). s 

As shown by the regression analysis, in most cases statistically significant for a 
particular character pair are linear and parabolic regressions. Out of 780 coefficients 
analysed only in 17 cases (2.2%) was an inconsistent result obtained. In 8 cases (1%) 
statistical significance of only the parabolic regression, and in 9 cases (1.2%) — of only 
the linear regression has been found. 

On the basis of the results from correlation and regression analyses it may be 
assumed that over the values analysed the relationship between most of the soil 
parameters studied can be reduced to a straight line. Linear regression for describing 
the interrelationships of the parameters studied is also suggested by the fact that if both 
regressions are significant, the respective parabole segments determining the type of 
correlation of the characters in question, over the value range analysed, are in fact 
rectilinear. One cannot, however, rule out the possibility that for a correlation and 
regression analysis of these parameters over a considerably wider range of values a 
curve of a higher-grade function would better represent their relationship. 

In the case of characters considered highly correlated, according to the assumption 
made, only the type of relationship found between magnesium and potassium, and 
manganese and soil ion exchange capacity (for the sub-Boreal variety) is approximated 
with a parabola. A parabola has also been adopted for defining the type of correlation 
between: magnesium and sodium, magnesium and total metal cation content, and 
magnesium and degree of saturation of the soil ion exchange capacity (for the 
sub-Boreal variety), because the values of correlation coefficients for parabolic 
correlation are higher (by about 40%) than linear correlation coefficients. 
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Fig. 3. Lines of regression of the characters 
a — total silt and clay and potassium content: I — y = 0.005 x + 0.121, II — y = 0.008 x + 0.037, II — 
y = 0.007 x + 0.127, IV — y = 0.006 x + 0.171, V — y = 0.014 x + 0.030, m.f. — mechanical fraction; 
b — calcium content and magnesium content: I — y = 0.105 x + 0.195, II — y = 0.094 x + 0.037, III — 
y = 0.079 x + 0.099, IV — y = 0.069 x + 0.589, V — y = 0.045 x + 0.540; c — content of potassium and 
content of sodium: | — y = 0.143 x + 0.023, Hl — y = 0.267 x + 0.019, III — y = 0.120 x + 0.025, IV — 
y = 0.185 x — 0.082, V — y = 0210 x — 0.013; d — soil pH and degree of soil ion exchange capacity 
saturation with metal cations: | — y = —0.032 x + 34.055, II — y = —0.032 x + 40.650, III — 
y = —0.046 x + 42.246, IV — y = —0.025 x + 37.103, V — y = —0.166 x + 59.092. All calculations for pH 
were based on specific acidity. | — sub-Boreal variety, LI — central-Polish variety, 111 — upland form of the 
Little-Polish variety, IV — piedmont form of the Little-Polish variety, V — Volhynian variety 
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Fig. 4. Lines of regression of the characters 
a — calcium content and total content of metal cations: I — y = 1.123 x + 0.285, II — y = 1.103 x + 0.207. 
Ml — y = 1.095 x + 0.392, IV — y = 1.085 x + 0.927, V — y = 1.071 x + 0.833; b — content of magnesium 
and total content of metal cations: 1 — y = —3.250 x? + 16.143 x — 2.216, II — y = 8.511 x + 1.685, 111 — 
y = 10931 x + 0.655, IV — y = 8.313 x + 0,619, V — y = 7.229 x + 3.820; c — content of potassium 
and total content of metal cations: I — y = 26.046 x + 0.499, III — y = 18.724 x + 1.184, IV — 
y = 18.838 x + 2.770, V — y = 29.596 x — 1.889; d — content of sodium and total content of metal 
cations: | — y = 123.036 x — 1.520, II — y = 110.650 x + 0.063, III — y = 107.980 x + 0.268, IV — 
y = 75.804 x + 5.727, V — y = 143.046 x — 0.272. I, II, III, IV, V — as in Figure 3 
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A regression analysis of the soil properties studied indicates that in the case of some 
characters the regression coefficients are similar in all the biotopes studied, and in the 
case of other characters — clearly different. Apart from this, for a particular character 
pair a more or less distinct regression coefficient similarity is seen only in some 
biotopes. The rate of increase in the value of y with unit changes in the value of x is 
graphically represented by the sloping of the regression line towards the axes of 
coordinates (Figs. 3, 4). 

The results that have been obtained suggest that the relationship between the 
parameters analysed better describes individual biotopes than a whole sample. This 
statement applies to most of the parameters taken into account in the analysis. 

Noteworthy are the results of the regression analysis describing the numeric ratios 
of total metal cation content and individual metal cation content (Fig. 4). In this case 
the correlation analysis and the description of the relationship with a regression 
equation have in fact made it possible to establish the proportions of the particular 
cations in the so-called total of bases (S). The course of the regression lines drawn for: 
total metal cation content and calcium content is almost the same in all the biotopes 
explored (Fig. 4 a). There are, therefore, grounds for finding an equation of regression 
of these characters for a whole sample: y = 1.103 x + 0.431. It makes it possible to 
estimate the total metal cation content from the content of calcium. 

The total content of metal cations can be estimated on the basis of calcium cation 
content also by evaluating the Ca-S~' expression. The adoption of this calculation is 
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Fig. 5. Lines of regression of the characters 
Total of magnesium, potassium and sodium cations, and the total content of metal cations: I — 
y = 8.261 x — 1.530, II — y = 7.672 x + 0.466, III — y = 8.264 x — 8.262, IV — y = 7.022 x — 0.256, V — 
y = 7.048 x + 0.195. I, II, III, IV, V — as in Figure 3 
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connected with the assumption that the proportion of calcium in the total metal cation 
content is constant for each S value. 

The percentages of calcium in the total metal cation content, calculated for the basic 
data, indicate that for low values of the total of bases calcium content ranges from 60 
to 70%, whereas for S values above 4 me: 100 g soil ^! it increases gradually, attaining 
a level of 80— 90%. Thus by estimating the total content of metal cations from the y 
= ax + b equation a more accurate estimation of the value of y is obtained than when 
the estimation is based on the Ca- S ^! expression. 

The proportions of the remainder of cations (i.e, Mg* *, K+, Na?) in the total 
content of metal cations, in all the biotopes studied, do not show a clear similarity like 
that found for calcium (Fig. 4 b, c, d). A high similarity is only found for the aggregate 
total of these cations in the five combinations analysed. For a whole sample it is 
expressed with the equation: y — 7.160 x — 0.241 (Fig. 5). 

The degree of correlation of calcium with the total content of the remainder of 
cations is high. The correlation coefficient value for these characters is 0.8, and the 
regression equation is of the format: y = 0.103 x + 0.429 (x — Ca^", y — Mg** + 
+K* + Na*). 

To sum up, the studies have shown that: 

(1) Most of the soil properties analysed in the present study are significantly 
positively correlated, but only some of them show a statistically significant interrela- 
tionship in all the five biotopes studied. 

(2) Over the range of values taken into account, the soil parameters are 
characterized by the simplest, i.e., linear type of relationship. Only in the case of very 
few characters is their correlation better described by a parabola. 

(3) Of the parameters taken into account in the analysis only the relationship 
between calcium cation content and the total metal cation content may characterize all 
the biotopes under study. A regression equation for these characters may also make it 
possible to estimate the total metal cation content from the content of calcium cations. 

(4) In the case of the remaining parameters the results from a regression analysis 
provide a better description of the soil environment of the particular oak-hornbeam 
forests studied than of an entire sample, that is to say, the habitat type analysed. 

(5) An analysis of stochastic relationships of the soil parameters in the diagram 
presented in the paper has made it possible to take into account soils that are different 
typologically and in respect of the advancement and intensity of the soil processes 
within the same soil type. The results thus relate to a fairly numerous group of soils 
which represent, broadly speaking, moist and mesotrophic habitats. 


4. SUMMARY 


The paper contains the results from correlation and regression analyses of 13 parameters determined in 
492 samples from the humus horizon (A,) of the soils of mixed oak-hornbeam forest communities of the 
Tilio-Carpinetum type. The data, providing a basis for a statistical analysis, characterize five biotope types 
representing oak-hornbeam forest habitats. 
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The correlation and regression analyses (by the least square method) have been carried out separately 
on the basis of numerical data for the individual biotopes corresponding to five typological units of mixed 
oak-hornbeam forests of the Tilio-Carpinetum type. Values of correlation coefficients of the characters 
analysed have been presented in Tables 1— 5. 

Most of the parameters taken into account in the analysis are significantly correlated with each other 
(Fig. 1). Characters considered highly correlated, according to the assumption in the paper, are in essence 
connected with the sorptive properties of the soil (Fig. 2). 

The results from the regression anaiysis indicate that over the range of values analysed the relationships 
between most soil parameters can be reduced to a straight line (Figs. 3— 5). 

The course of the regression lines describing the numeric ratios of: total content of metal cations and 
calcium content provides a basis for finding an equation of regression of these characters for a whole sample 
(Fig. 4a). This equation: y = 1.103 x + 0.431 can be used for the estimation of the so-called total of bases from 
the content of calcium cations. 

In the case of the remaining characters the results from the analysis of regression provide a better 
description of the particular biotopes than of a whole sample, i.e., the habitat type analysed. 


5. POLISH SUMMARY 


W pracy zawarto wyniki analizy korelacji i regresji 13 parametrów oznaczonych w 492 próbach z 
poziomu akumulacyjno-próchnicznego (A,) gleb zbiorowisk gradowych typu Tilio-Carpinetum. Dane 
bedace podstawa analizy statystycznej charakteryzuja 5 typów biotopów reprezentujacych siedliska lasów 
gradowych. 

Analize korelacji i regresji (metoda najmniejszych kwadratów) wykonano oddzielnie na podstawie 
danych liczbowych dla poszczególnych biotopów odpowiadajacych 5 jednostkom typologicznym gradów 
typu Tilio-Carpinetum. Wartosci wspólczynników korelacji analizowanych cech zawieraja tab. 1 — 5. 

Większość uwzględnionych w analizie parametrów jest ze sobą istotnie skorelowana (rys. 1). Cechy. 
które zgodnie z przyjętym w pracy założeniem uznano za silnie sprzężone, są w zasadzie związane z 
właściwościami sorpcyjnymi gleby (rys. 2). - 

Wyniki analizy regresji wskazują, że w zakresie analizowanych wartości, zależność między większością 
parametrów glebowych można sprowadzić do linii prostej (rys. 3— 5). 

Przebieg linii regresji opisujących stosunki ilościowe cech: ogólna zawartość kationów metali i zawartość 
wapnia daje podstawy do wyznaczenia równania regresji tych cech dla calej próby (rys. 4 a). Równanie to, 
y = 1,103 x + 0.431, może służyć do szacowania wartości tzw. sumy zasad na podstawie zawartości 
kationów wapnia. 

W przypadku pozostałych cech wyniki analizy regresji lepiej charakteryzują poszczególne biotopy niż 
całą próbę. tzn. analizowany typ siedliska. 
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